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TITLE OF THE INVENTION 

A Vehicle Controller for Controlling an Air-fuel Ratio 

BACKGROUND OF THE INVENTION 
5 Technical Field 

The invention relates to a controller for controlling an air-fuel ratio 
based on an output of an exhaust gas sensor disposed in an exhaust system 
of an internal-combustion engine. 

10 Description of the Related Art 

A catalyst converter is provided in an exhaust system of an internal 
combustion engine of a vehicle. When the air-fuel ratio of air-fuel mixture 
introduced into the engine is lean, the catalyst converter oxidizes HC and 
CO with excessive oxygen included in the exhaust gas. When the air-fuel 

15 ratio is rich, the catalyst converter reduces Nox with HC and CO. When the 
air-fuel ratio is in the stoichiometric air-fuel ratio region, HC, CO and Nox 
are simultaneously and effectively purified. 

An exhaust gas sensor is provided downstream of the catalyst 
converter. The exhaust gas sensor detects the concentration of oxygen 

20 included in the gas that is discharged into the exhaust system. Feedback 
control for the air-fuel ratio of the engine is performed based on the output 
of the exhaust gas sensor. 

As an example of the feedback control for the air-fuel ratio, 
Japanese Patent Application Unexamined Publication No. Hll- 153051 

25 proposes response assignment control in which a switching function is 
defined. This control converges the output of the exhaust gas sensor to a 
target value by converging the value of the switching function to zero. A 
controlled variable (a target air-fuel ratio) for converging the output of the 
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exhaust gas sensor to the target value is calculated. The output of the 
exhaust gas sensor and the output of the air-fuel ratio (LAF) sensor, which 
is disposed upstream of the catalyst converter for detecting the air-fuel 
ratio, are used for the calculation of the controlled variable. A fuel amount 
to be supplied to the engine is controlled according to the calculated 
controlled variable. 

Recently, there is a trend to enhance a response of the exhaust gas 
sensor so as to stabilize the accuracy of detecting deterioration of the 
catalyst and to decrease the amount of discharged NOx. When a response of 
the exhaust gas sensor is enhanced, high-frequency components, which are 
called chemical noise, may be introduced into the output of the exhaust gas 
sensor. Such chemical noise may cause variations in the target air-fuel 
ratio because the target air-fuel ratio is calculated based on the output of 
the exhaust gas sensor. Such variations in the target air-fuel ratio may 
cause a large variation in the actual air-fuel ratio, which reduces the 
purification rate of the catalyst. 

The air-fuel ratio is sometimes made rich so as to protect the engine 
and the catalyst. Such enrichment of the air-fuel ratio increases the 
amount of discharged CO. In order to suppress the discharge of CO, it is 
preferable to perform the air-fuel ratio control in a form of closed loop. On 
the other hand, the air-fuel ratio control for making the air fuel ratio lean 
may be performed so as to improve the fuel efficiency. In such a state in 
which the air-fuel ratio is made lean, it is preferable to perform the air-fuel 
ratio control in the form of closed loop. In order to stably perform the 
closed-loop air-fuel ratio control, there is a trend to expand a detection 
range of the air-fuel ratio (LAF) sensor. 

There is a limitation in the performance of an A/D converter that 
converts an analog signal from the air-fuel ratio sensor into a digital signal. 
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When a detection range of the air-fuel ratio sensor is expanded, such 
limitation of the performance of the A/D converter reduces the resolution of 
the air-fuel ratio detected by the air-fuel ratio sensor. Such resolution 
reduction may reduce the capability to make the actual air-fuel ratio follow 
the target air-fuel ratio in the air-fuel ratio control, which reduces the 
purification rate of the catalyst. Such resolution reduction may also reduce 
the accuracy of identifying a model parameter for the air-fuel ratio control 
because the model parameter is identified based on the actual air-fuel ratio. 
The reduction of the accuracy of identifying a model parameter may also 
reduce the purification rate of the catalyst. 

Therefore, there is a need for an apparatus and a method capable of 
removing chemical noise from the output of the exhaust gas sensor when a 
response of the exhaust gas sensor is enhanced. There is also a need for an 
apparatus and a method capable of compensating the shortage of resolution 
of the air-fuel ratio sensor when a detection range of the air-fuel ratio 
sensor is expanded. 

SUMMARY OF THE INVENTION 

According to a first embodiment of the invention, a controller for 
controlling an air-fuel ratio of an internal-combustion engine comprises a 
first exhaust gas sensor for detecting oxygen concentration of the exhaust 
gas, a first decimation filter connected to the first exhaust gas sensor, and a 
control unit connected to the first decimation filter. The control unit 
determines a manipulated variable for manipulating the air-fuel ratio so 
that an output value from the first decimation filter converges to a target 
value. The first decimation filter further comprises a first oversampler, a 
first low-pass filter, and a first downsampler. The first oversampler 
oversamples the output of the first exhaust gas sensor in a shorter cycle 
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than a cycle that is used for determining the manipulated variable. The 
first low-pass filter smoothes the oversampled value. The first 
downsampler re-samples the smoothed value in the cycle that is used for 
determining the manipulated variable. Thus, the first decimation filter 
outputs the re-sampled value. 

As described above, when a response of the first exhaust gas sensor 
is enhanced, chemical noise may appear in the output of the first exhaust 
gas sensor. The first decimation filter can remove such chemical noise. The 
air-fuel ratio control based on the output from the first decimation filter 
prevents the purification rate of the catalyst from deteriorating. 

According to a second embodiment of the invention, a controller for 
controlling an air-fuel ratio of an internal-combustion engine comprises a 
first exhaust gas sensor provided downstream of a catalyst converter, a 
second exhaust gas sensor provided upstream of the catalyst converter, a 
second decimation filter connected to the second exhaust gas sensor, and a 
control unit connected to the second decimation filter. The first exhaust gas 
sensor detects oxygen concentration of the exhaust gas. The second exhaust 
gas sensor detects an air-fuel ratio of the exhaust gas. The control unit uses 
an output value from the second decimation filter to determine a 
manipulated variable for manipulating the air-fuel ratio so that an output 
value from the first exhaust gas sensor converges to a target value. The 
second decimation filter further comprises a second oversampler, a second 
low-pass filter, and a second downsampler. The second oversampler 
oversamples the output of the second exhaust gas sensor in a shorter cycle 
than a cycle that is used for determining the manipulated variable. The 
second low-pass filter smoothes the oversampled value. The second 
downsampler re-samples the smoothed value in the cycle that is used for 
determining the manipulated variable. Thus, the second decimation filter 
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outputs the re-sampled value. 

As described above, when a detection range of the second exhaust 
gas sensor is expanded, the resolution of the air-fuel ratio detected by the 
second exhaust gas sensor may be reduced. The second decimation filter 
can compensate the shortage of resolution of the second exhaust gas sensor. 
Specifically, the second decimation filter estimates detection values below 
the resolution limit of the second exhaust gas sensor. The air-fuel ratio 
control based on the output from the second decimation filter prevents the 
purification rate of the catalyst from deteriorating. 

According to one embodiment of the invention, the manipulated 
variable is determined by response assignment control. The response 
assignment control can stably and quickly cause the output of the first 
exhaust gas sensor to converge to a target value. 

According to another embodiment of the invention, the manipulated 
variable is determined by performing control that uses one of A£ 
(delta-sigma) modulation algorithm, A (delta) modulation algorithm and £A 
(sigma-delta) modulation algorithm. The control using such an algorithm 
can stably cause the output of the first exhaust gas sensor to converge to a 
target value even when a delay in the response of an object to be controlled 
by the air-fuel ratio control is large. 

According to one embodiment of the invention, an object to be 
controlled is an exhaust system. The exhaust system extends from the 
second exhaust gas sensor through the catalyst converter to the first 
exhaust gas sensor. In the first embodiment described above, a dead time in 
the exhaust system is determined based on the output value from the first 
decimation filter. An estimated value for the output of the first exhaust gas 
sensor is calculated so that the dead time in the exhaust system is 
compensated. The estimated value is used to determine the manipulated 
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variable for manipulating the air-fuel ratio. Since the estimated value for 
the output of the first exhaust gas sensor is determined considering the 
dead time in the exhaust system, the manipulated variable enables the 
output value from the first exhaust gas sensor to stably converge to a target 
value. In the second embodiment described above, the dead time in the 
exhaust system is determined based on the output value from the second 
decimation filter. 

According to yet another embodiment of the invention, the object of 
the air-fuel ratio control further includes an air-fuel ratio manipulating 
system. The air-fuel ratio manipulating system extends from the control 
unit for determining the manipulated variable through the engine to the 
second exhaust gas sensor. In the first embodiment described above, a dead 
time in the air-fuel ratio manipulating system is determined based on the 
output value from the first decimation filter. An estimated value for the 
output of the first exhaust gas sensor is calculated so that the dead time in 
the exhaust system and the dead time in the air-fuel ratio manipulating 
system are compensated. The estimated value is used to determine the 
manipulated variable for manipulating the air-fuel ratio. Since the 
estimated value for the output of the first exhaust gas sensor is determined 
considering the dead time both in the exhaust system and in the air-fuel 
ratio manipulating system, the manipulated variable enables the output 
value from the first exhaust gas sensor to stably converge to a target value. 
In the second embodiment described above, the dead time in the air-fuel 
ratio manipulating system is determined based on the output value from 
the second decimation filter. 

According to yet another embodiment of the invention, the control 
unit calculates a parameter that is used for determining the manipulated 
variable. The parameter acts to adapt the air-fuel ratio manipulation to 
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state changes of the exhaust system. In the first embodiment described 
above, the parameter is calculated based on the output value from the first 
decimation filter. Since the output from the first decimation filter does not 
include chemical noise, the parameter is calculated with a better accuracy. 
In the second embodiment described above, the parameter is calculated 
based on the output value from the second decimation filter. Since the 
second decimation filter provides detection values below the resolution 
limit of the second exhaust gas sensor, the parameter is calculated with a 
better accuracy. 

According to yet another embodiment of the invention, a cut-off 
frequency of the first and second low-pass filters of the first and second 
decimation filters is set to a higher frequency than a frequency that is used 
for detecting a failure of the catalyst. Thus, the air-fuel ratio control can be 
performed without reducing the accuracy of detecting a failure of the 
catalyst. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a schematic view of an internal combustion engine and its 
controller according to one embodiment of the present invention. 
Fig. 2 is a view of layout of a catalyst converter and an exhaust gas sensor 
according to one embodiment of the present invention. 

Figure 3 shows an outline of air-fuel ratio control according to one 
embodiment of the present invention. 

Figure 4 is a block diagram of air-fuel ratio control according to a first 
embodiment of the present invention. 

Figure 5 is a detailed functional block diagram of a controller according to 
one embodiment of the present invention. 

Figure 6 schematically shows a switching line for response assignment 
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control according to one embodiment of the present invention. 

Figure 7 shows response characteristics of response assignment control 

according to one embodiment of the present invention. 

Figure 8 is a detailed functional block diagram of a first decimation filter 
according to one embodiment of the present invention. 

Figure 9 shows a frequency response of an output of an exhaust gas sensor 
that varies according to the degree of deterioration of catalyst. 
Figure 10 shows a frequency response of a filtered output of an exhaust gas 
sensor that varies according to the degree of deterioration of catalyst. 
Figure 11 shows low-pass filter characteristics of a first decimation filter 
according to one embodiment of the present invention. 

Figure 12 shows a filtered output of an exhaust gas sensor (a) when a first 
decimation filter is used, and (b)(c) when another filter is used, according to 
one embodiment of the present invention. 

Figure 13 shows low-pass filter characteristics of a filter that is used in the 
case (b) of Figure 12. 

Figure 14 shows low-pass filter characteristics of a filter that is used in the 
case (c) of Figure 12. 

Figure 15 shows an appearance of chemical noise in the exhaust gas sensor 
output, and variations in the target air-fuel ratio error kcmd according to 
conventional air-fuel ratio control. 

Figure 16 shows the output of the first decimation filter Vo2_df and the 
target air-fuel ratio error kcmd according to one embodiment of the present 
invention. 

Figure 17 is a block diagram of air-fuel ratio control according to a second 
embodiment of the present invention. 

Figure 18 is a detailed functional block diagram of a second decimation 
filter according to one embodiment of the present invention. 

8 
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Figure 19 shows low-pass filter characteristics of a second decimation filter 
according to one embodiment of the present invention. 

Figure 20 shows a shortage of resolution of an air-fuel ratio (LAF) sensor, 
the target air-fuel ratio error kcmd, and an identified parameter "bl" 
according to conventional air-fuel ratio control. 

Figure 21 shows the output of a second decimation filter KACT_df and the 
target air-fuel ratio error kcmd according to one embodiment of the present 
invention. 

Figure 22 is a block diagram of air-fuel ratio control according to third 
embodiment of the present invention. 

Figure 23 shows a detailed functional block diagram of a controller 
according to one embodiment of the present invention. 

Figure 24 shows a detailed functional block diagram of a DMS controller 
according to one embodiment of the present invention. 

Figure 25 is a flowchart of a main routine for adaptive air-fuel ratio control 
according to one embodiment of the present invention. 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 
Structure of internal-combustion engine and control apparatus 

Preferred embodiments of the present invention will be described 
referring to the attached drawings. Figure 1 is a block diagram showing a 
controller of an internal-combustion engine (hereinafter referred to as an 
engine) in accordance with one embodiment of the invention. 

An electronic control unit (hereinafter referred to as an ECU) 5 
comprises an input interface 5a for receiving data sent from each part of 
the engine 1, a CPU 5b for carrying out operations for controlling each part 
of the engine 1, a storage device 5c including a read only memory (ROM) 
and a random access memory (RAM), and an output interface 5d for 

9 
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sending control signals to each part of the engine 1. Programs and various 
data for controlling each part of the vehicle are stored in the ROM. A 
program for controlling an air-fuel ratio according to the invention, data 
and tables used for operations of the program are stored in the ROM. The 
5 ROM may be a rewritable ROM such as an EEPROM. The RAM provides 
work areas for operations by the CPU 5a, in which data sent from each part 
of the engine 1 as well as control signals to be sent out to each part of the 
engine 1 are temporarily stored. 

The engine 1 is, for example, an engine equipped with four cylinders. 
10 An intake manifold 2 is connected to the engine 1. A throttle valve 3 is 
disposed upstream of the intake manifold 2. A throttle valve opening ( d TH) 
sensor 4, which is connected to the throttle valve 3, outputs an electric 
signal corresponding to an opening angle of the throttle valve 3 and sends 
it to the ECU 5. 

15 A bypass passage 21 for bypassing the throttle valve 3 is provided in 

the intake manifold 2. A bypass valve 22 for controlling the amount of air to 
be supplied into the engine 1 is provided in the bypass passage 21. The 
bypass valve 22 is driven in accordance with a control signal from the ECU 
5. 

20 A fuel injection valve 6 is provided for each cylinder at an 

intermediate point in the intake manifold 2 between the engine 1 and the 
throttle valve 3. The fuel injection valve 6 is connected to a fuel pump (not 
shown) to receive fuel supplied from a fuel tank (not shown). The fuel 
injection valve 6 is driven in accordance with a control signal from the ECU 

25 5. 

An intake manifold pressure (Pb) sensor 8 and an outside air 
temperature (Ta) sensor 9 are mounted in the intake manifold 2 
downstream of the throttle valve 3. The detected intake manifold pressure 
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Pb and outside air temperature Ta are sent to the ECU 5. 

An engine water temperature (TW) sensor 10 is attached to the 
cylinder peripheral wall, which is filled with cooling water, of the cylinder 
block of the engine 1. The temperature of the engine cooling water detected 
by the TW sensor is sent to the ECU 5. 

A rotational speed (Ne) sensor 13 is attached to the periphery of the 
camshaft or the periphery of the crankshaft (not shown) of the engine 1, 
and outputs a CRK signal pulse at a predetermined crank angle cycle (for 
example, a cycle of 30 degrees) that is shorter than a TDC signal pulse 
cycle issued at a crank angle cycle associated with a TDC position of the 
piston. CRK pulses are counted by the ECU 5 to determine the rotational 
speed Ne of the engine 1. 

An exhaust manifold 14 is connected to the engine 1. The engine 1 
discharges exhaust gas through the exhaust manifold 14. A catalyst 
converter 15 removes deleterious substances such as HC, CO, and Nox 
included in exhaust gas flowing through the exhaust manifold 14. The 
catalyst converter 15 comprises two catalysts, an upstream catalyst and a 
downstream catalyst. 

A full range air-fuel ratio (LAF) sensor 16 is provided upstream of 
the catalyst converter 15. The LAF sensor 16 linearly detects the 
concentration of oxygen included in exhaust gas over a wide air-fuel ratio 
zone, from the rich zone where the air/fuel ratio is richer than the 
stoichiometric air/fuel ratio to an extremely lean zone. The detected oxygen 
concentration is sent to the ECU 5. 

An 02 (exhaust gas) sensor 17 is provided between the upstream 
catalyst and the downstream catalyst. The 02 sensor 17 is a binary -type of 
exhaust gas concentration sensor. The 02 sensor outputs a high level signal 
when the air-fuel ratio is richer than the stoichiometric air-fuel ratio, and 

11 
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outputs a low level signal when the air-fuel ratio is leaner than the 
stoichiometric air-fuel ratio. The electric signal is sent to the ECU 5. 

The 02 sensor 17may be referred to as a first exhaust gas sensor. 
The LAF sensor 16 may be referred to as a second exhaust gas sensor. 
5 Signals sent to the ECU 5 are passed to the input circuit 5a. The 

input interface 5a converts analog signal values into digital signal values. 
The CPU 5b processes the resulting digital signals, performs operations in 
accordance with the programs stored in the ROM, and creates control 
signals. The output interface 5d sends these control signals to actuators for 
10 the bypass valve 22, fuel injection valve 6 and other mechanical 
components. 

Figure 2 shows a structure of the catalyst converter 15. Exhaust gas 
introduced into the exhaust manifold 14 passes through the upstream 
catalyst 25 and then through the downstream catalyst 26. It is known that 

15 it is easier to maintain the purification rate of Nox at an optimal level by 
air-fuel ratio control based on the output of an 02 sensor provided between 
the upstream and downstream catalysts, compared with air-fuel ratio 
control based on the output of an 02 sensor provided downstream of the 
downstream catalyst. Therefore, in the embodiment of the invention 

20 described hereafter, the 02 sensor 17 is provided between the upstream 
and downstream catalysts. The 02 sensor 17 detects the concentration of 
oxygen included in the exhaust gas after the passage through the upstream 
catalyst 25. 

Alternatively, the 02 sensor may be disposed downstream of the 
25 downstream catalyst 26. If the catalyst converter 15 is implemented with a 
single catalyst, the 02 sensor is disposed downstream of the catalyst 
converter 15. 

Figure 3 shows purification behavior of the upstream catalyst and 
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the downstream catalyst. A window 27 indicates an air-fuel ratio region in 
which CO, HC and Nox are optimally purified. Since oxygen included in 
exhaust gas is consumed by the purification in the upstream catalyst 25, 
the exhaust gas supplied to the downstream catalyst 26 exhibits a 
5 reduction atmosphere (i.e., a rich state) as shown by a window 28. In such a 
reduction atmosphere, Nox is further purified. Thus, the cleaned exhaust 
gas is discharged. 

In order to optimally maintain the purification performance of the 
catalyst converter 15, adaptive control of the air-fuel ratio according to the 
10 invention causes the output of the 02 sensor 17 to converge to a target 
value so that the air-fuel ratio is within the window 27. 

A reference number 29 shows an allowable range that defines a 
limitation of a variable manipulated by the adaptive air-fuel ratio control, 
which will be described in detail later. 

15 

Adaptive air-fuel ratio control in accordance with a first embodiment 

Figure 4 shows a block diagram of adaptive air-fuel ratio control in 

accordance with a first embodiment of the invention. The LAF sensor 16 

detects an air-fuel ratio KACT of the exhaust gas supplied to the upstream 
20 catalyst 25. The 02 sensor 17 outputs a voltage Vo2/OUT that indicates the 

oxygen concentration of the exhaust gas after the purification by the 

upstream catalyst 25. 

The output Vo2/OUT from the 02 sensor 17 is delivered to a first 

decimation filter 36. The first decimation filter 36 oversamples the output 
25 V62/OUT of the 02 sensor 17, performs a low-pass filtering process on the 

oversampled value, and then downsamples the filtered value. The output 

from the first decimation filter 36 is represented by a sample value Vo2_df. 

The sample value Vo2_df is compared with a target value Vo2/TARGET. An 
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error Vo2 between the sample value Vo2_df and the target value 
V62/TARGET is supplied to a controller 31. 

An object (i.e., plant) to be controlled by the adaptive air-fuel ratio 
control is an exhaust system 19 extending from the LAF sensor 16 through 
5 the upstream catalyst 25 to the 02 sensor 17. The controller 31 determines 
a target air-fuel ratio error "kcmd" based on the error Vo2. The target 
air-fuel ratio error kcmd is added to a base value FLAF/BASE to determine 
a target air-fuel ratio KCMD. A fuel injection amount is corrected 
accordance with the target air-fuel ratio KCMD and is supplied to the 
10 engine 1. After the fuel injection, the output Vo2/OUT of the 02 sensor 17 is 
detected again. 

Thus, the controller 31 performs a feedback control to determine the 
target air-fuel ratio KCMD so that the error Vo2 converges to zero. The 
exhaust system 19, which is the controlled object, can be modeled as shown 
15 by the equation (l) in which Vo2/OUT is defined as a control output and the 
output KACT of the LAF sensor is defined as a control input. The exhaust 
system 19 is modeled as a discrete-time system. Such modeling can make 
the air-fuel ratio control algorithm simple and suitable for computer 
processing. 

20 

Vo2 ? (k + 1) = al ■ Vo2'(k) + a2 ■ Vo2'(k - 1) + bl • kact(k - d3) 
where Vo2' (k) = Vo2 / OUT(k) - Vo2 / TARGET 
kact(k) = KACT(k) - FLAF/BASE 

(1) 

"k" is an identifier for identifying a control cycle. A sensor output 
error Vo2' indicates an error between the 02 sensor output V62/OUT and 
the target value Vo2/TARGET. An actual air-fuel ratio error "kact" indicates 
25 an error between the LAF sensor output KACT and the base value 
FLAF/BASE. The base value FLAF/BASE is set to be a central value for the 
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target air-fuel ratio KCMD. For example, the base value is set to a value 
indicative of stoichiometry (that is, FLAF/BASE=l). The base value 
FLAF/BASE may be a constant value, or may be established according to 
the operating state of the engine. 
5 "d3" indicates a dead time in the exhaust system 19. The dead time 

d3 is a time required for the air-fuel ratio detected by the LAF sensor 16 to 
be reflected in the output of the 02 sensor 17. "al", "a2" and "bl" are model 
parameters, which are generated by an identifier. The identifier will be 
described later. 

10 Relation between the 02 sensor output Vo2/OUT and the output 

Vo2_df of the first decimation filter is expressed as shown in the equation 
(2). 

Vo2 / OUT(k) = Vo2 _ df (k + d5) 

(2) 



"d5" indicates a dead time in the first decimation filter 36. The dead 

15 time d5 is a time required for the 02 sensor output Vo2/OUT to be 

oversampled, filtered using a low-pass filter and then downsampled. The 

dead time d5 is, for example, one control cycle (that is, d5=l). 

A system including the exhaust system 19 and the first decimation 

filter 36 is determined based on the equations (l) and (2), as shown in the 

20 equation (3). 

Vo2(k + 1) = al • Vo2(k) + a2 • Vo2(k - 1) + bl • kact(k - d3 - d5) 
= al • Vo2(k) + a2 • Vo2(k - 1) + bl • kact(k - dl) 
where Vo2(k) = Vo2 _ df (k) - Vo2 / TARGET 
kact(k) = KACT(k) - FLAF / BASE 
dl = d3 + d5 

(3) 

Thus, incorporation of the first decimation filter 36 increases the 

15 
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dead time in the exhaust system. 

On the other hand, an air-fuel ratio manipulating system 18 
extending from the ECU 5 through the engine 1 to the LAF sensor 16 can 
be modeled as shown by the equation (4). 

kact(k) = kcmd(k-d2) 

where kact(k) = KACT(k) - FLAF / BASE 
kcmd(k) = KCMD(k) - FLAF /BASE 

(4) 

The target air-fuel ratio error "kcmd" indicates an error between the 
target air-fuel ratio KCMD and the base value FLAF/BASE. "d2" indicates 
a dead time in the air-fuel ratio manipulating system 18. The dead time d2 
is a time required for the calculated target air-fuel ratio KCMD to be 
reflected in the output KACT of the LAF sensor 16. 

The air-fuel ratio manipulating system 18 may be included in the 
object to be controlled by the adaptive air-fuel ratio control. In this case, 
the model equation is expressed based on the equations (3) and (4), as 
shown by the equation (5). A dead time "d" is a total dead time in a system 
comprising the air-fuel ratio manipulating system 18, the exhaust system 
19, and the first decimation filter 36. Incorporation of the first decimation 
filter 36 increases the dead time. 

Vo2(k + 1) = al ■ Vo2(k) + a2 • Vo2(k - 1) + bl - kcmd(k - dl - d2) 
= al • Vo2(k) + a2 ■ Vo2(k - 1) + bl • kcmd(k - d) 
where d = dl + d2 = d3 + d5 + d2 

(5) 

Figure 5 shows a more detailed block diagram of the controller 31 
shown in Figure 4. The controller 31 comprises an identifier 32, an 
estimator 33, a sliding mode controller 34, and a limiter 35. 

The identifier 32 identifies the model parameters al, a2 and bl in 

16 
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the equation (3) so that a modeling error is removed. An identification 
method performed by the identifier 32 will be described. 

The identifier 32 uses model parameters al(k-l), a2(k-l) and 
bl(k-l) that have been calculated in the previous control cycle to 
5 determine a sensor output error V62(k) for the current cycle in accordance 
with the equation (6). 

V62(k) = al(k - 1) • Vo2(k - 1) 
+ a2(k-l)-Vo2(k-2) 
+ bl(k-l)-kact(k-dl-l) 

(6) 

The equation (7) indicates an error id/e(k) between the sensor 
output error V62(k) that is calculated in accordance with the equation (6) 
10 and a sensor output error Vo2(k) that is actually detected in the current 
control cycle. 

id / e(k) = Vo2(k) - V62(k) (7) 



15 The identifier 32 calculates al(k), a2(k) and bl(k) for the current 

cycle so that the error id/e(k) is minimized. Here, a vector 9 is defined as 
shown in the equation (8). 



20 



0 T (k) = [al(k) a2(k) bl(k)] (8) 



The identifier 32 determines al(k), a2(k) and bl(k) in accordance 
with the equation (9). As shown by the equation (9), al(k), a2(k) and bl(k) 
for the current control cycle are calculated by changing al(k), a2(k) and 
bl(k) calculated in the previous control cycle by an amount proportional to 
25 the error id/e(k). 
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0(k) = 0(k - 1) + K9(k) • id / e(k) 



where K0(k) = 



P(k-l)4(k) 



l + ^ T (k)P(k-l)§(k) 
^ T (k) = [Vo2(k-l) Vo2(k-2) kact(k-dl-l)] 
1 - X2(k)P(k-l)4(k)^ T (k) 



P(k) = 



[I- 



M(k) L ~ M(k) + X2(k)^ T (k)P(k - l)5(k) 
0 < M < 1 0<X2<2 I: unit matrix 

(9) 



F(k-l) 



10 



In order to compensate the dead time "dl" of the exhaust system 19 
and the dead time "d2" of the air-fuel ratio manipulating system, the 
estimator 33 estimates a sensor output error Vo2 after the dead time d 
(=dl+d2) based on the model equation (5). Specifically, the estimated value 



Vo2(k + d) for the sensor output error Vo2(k+d) after the dead time "d" is 
determined, as shown by the equation (10). Coefficients al, a2 and (3 are 
calculated using model parameters determined by the identifier 32. Past 
time-series data kcmd(k-j) (wherein, j=l, 2, ... d) of the air-fuel ratio error 
includes air-fuel ratio errors obtained during a period of the dead time "d." 



Vo2(k + d) = al • Vo2(k) + a2 • Vo2(k - 1) + £ p"j • kcmd(k - j) 

j=i 

where al = first - row, first - column element of A d 

a2 = first - row, sec ond - column element of A d 
Pj = first row elements of A J_1 • B 
al a2" 



A = 



B = 



1 0 
bl~ 



(10) 



15 Past values kcmd(k-d2), kcmd(k-d2 l), ... kcmd(k-d) for the target 

air-fuel ratio error "kcmd" before the dead time d2 can be replaced with 
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actual air-fuel ratio errors kact(k), kact(k-l), ... kact(k-d+d2) by using the 
equation (2). As a result, the equation (ll) is derived. 



Vo2(k + d) = ctl • Vo2(k) + a2 • Vo2(k - 1) 

d2-l d-d2 

+ ]T (3j - kcmd(k - j) + £ pi + d2 • kact(k - i) 

j=l i=0 

= al • Vo2(k) + a2 • Vo2(k - 1) 

d2-l dl 

+ £ pj • kcmd(k - j) + 2 Pi + d2 • kact(k - i) 

j=l i=0 

(11) 

5 

The sliding mode controller 34 establishes a switching function a so 
as to perform the sliding mode control, as shown in the equation (12). 

a(k) = s ■ Vo2(k - 1) + Vo2(k) (12) 

10 

Vo2(k-l) indicates the sensor output error detected in the previous 
cycle as described above. Vo2(k) indicates the sensor output error detected 
in the current cycle, "s" is a setting parameter of the switching function a, 
and is established to satisfy -Ks<l. 
15 The equation in the case of a(k)=0 is called an equivalent input 

system, which specifies the convergence characteristics of the sensor output 
error Vo2, or a controlled variable. Assuming <r(k)=0, the equation (12) is 
transformed to the equation (13). 



20 Vo2(k) = -s - Vo2(k - 1) (1 3) 



Now, characteristics of the switching function a will be described 
with reference to Figure 6 and the equation (13). In Figure 6, the equation 
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(13) is shown as a line 41 on a phase plane with Vo2(k-l) being the 
horizontal axis and Vo2(k) being the vertical axis. The line 41 is referred to 
as a switching line. It is assumed that the initial value of a state variable 
(Vo2(k-l), Vo2(k)) that is a combination of Vo2(k*l) and Vo2(k) is shown by a 
5 point 42. The sliding mode control operates to place the state variable 
shown by the point 42 on the line 41 and then restrain it on the line 41. 
According to the sliding mode control, since the state variable is held on the 
switching line 41, the state variable can highly stably converge to the 
origin 0 of the phase plane without being affected by disturbances or the 

10 like. In other words, by confining the state variable (Vo2(k-l), Vo2(k)) on 
such a stable system having no input as shown by the equation (13), the 
sensor output error Vo2 can converge to zero robustly against disturbances 
and modeling errors. 

The switching function setting parameter "s" is a parameter which 

15 can be variably selected. Reduction (convergence) characteristics of the 
sensor output error Vo2 can be specified by the setting parameter "s." 

Figure 7 shows one example of response assignment characteristics 
of the sliding mode control. A line 43 shows a case in which the value of the 
setting parameter is "1." A curve 44 shows a case in which the value of the 

20 setting parameter is "0.8." A curve 45 shows a case in which the value of 
the setting parameter is "0.5." As seen from the figure, the rate of 
convergence of the sensor output error Vo2 changes according to the value 
of the setting parameter "s." It is seen that the convergence rate becomes 
faster as the absolute value of "s" becomes smaller. 

25 Three control inputs are determined to cause the value of the 

switching function a to converge to zero. That is, a control input Ueq for 
confining the state variable on the switching line, a control input Urch for 
placing the state variable on the switching line, and a control input Uadp 
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for placing the state variable on the switching line while suppressing 
modeling errors and disturbances. The three control inputs Ueq, Urch and 
Uadp are summed to determine a demand error Usl. The demand error Usl 
is used to calculate the air-fuel ratio error kcmd. 
5 The equivalent control input Ueq needs to satisfy the equation (14) 

because it is an input for restraining the state variable onto the switching 
line. 

a(k + l) = o(k) (14) 

10 

The equivalent control input Ueq that satisfies a(k+l)=o(k) is 
determined from the equations (4), (5) and (12), as shown by the equation 
(15). 

Ueq(k) = -— [((al - 1) + s) • Vo2(k + d) + (a2 - s) • Vo2(k + d - 1)] 
15 bl 

(15) 

The reaching law input Urch has a value that depends on the value 
of the switching function a. The reaching law Urch is determined in 
accordance with the equation (16). In the embodiment, the reaching law 
20 input Urch has a value proportional to the value of the switching function a. 
Krch indicates a feedback gain of the reaching law, which is predetermined 
with, for example, simulation in which the stability and quick response of 
convergence of the value of the switching function to zero (a=0) are taken 
into consideration. 

25 
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Urch(k) = -— • Krch ■ a(k + d) (16) 
bl 

The adaptive law input Uadp has a value that depends on an 
integrated value of the switching function a. The adaptive law input Uadp 
5 is determined in accordance with the equation (17). In the embodiment, the 
adaptive law input Uadp has a value proportional to the integrated value of 
the switching function a. Kadp indicates a feedback gain of the adaptive 
law, which is predetermined with, for example, simulation in which the 
stability and quick response of convergence of the value of the switching 
10 function to zero (<j=0) are taken into consideration. AT indicates the period 
of a control cycle. 

Uadp (k) = -^-- Kadp -£(a(i) -AT) (17) 
bl 7^ 

15 Since the sensor output errors Vo2(k+d) and Vo2(k+d-l), and the 

value a(k+d) of the switching function include the dead time "d", these 
values can not be directly obtained. Therefore, the equivalent control input 
Ueq is determined using an estimated errors Vo2(k + d) and Vo2(k-f d-1) 
generated by the estimator 33. 

20 

Ueq(k) = [((al - 1) + s) • Vo2(k + d) + (a2 - s) ■ Vo2(k + d - 1)] 
bl 

(18) 

A switching function a is determined using the estimated errors 
generated by the estimator 33, as shown in the equation (19). 
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a = s • Vo2(k - 1) + Vo2(k) (1 9) 

The switching function a is used to determine the reaching law 
5 input Urch and the adaptive law input Uadp. 

Urch(k) = - — • Krch • a(k + d) (20) 
bl 

i k+d _ 

Uadp(k) = -± ■ Kadp- £(a(i) ■ AT) (21) 
bl i= o 

10 

As shown by the equation (22), the equivalent control input Ueq, the 
reaching law input Urch and the adaptive law input Uadp are added to 
determine a demand error Usl. 

1 5 Usl(k) = Ueq(k) + Urch(k) + Uadp(k) (22) 

The limiter 35 performs a limiting process for the demand eror Usl 
to determine the air-fuel ratio error kcmd. More specifically, if the demand 
error Usl is within an allowable range, the limiter 35 sets the air-fuel ratio 

20 error kcmd to the value of the demand error Usl. If the demand error Usl 
deviates from the allowable range, the limiter 35 sets the air-fuel ratio 
error kcmd to an upper or lower limit value of the allowable range. 

As shown by reference number 29 in Figure 3, the allowable range 
used by the limiter 35 is set to a range whose center is almost located in the 

25 window 27 and whose width is wider than that of the window 27. The 
allowable range is actively established in accordance with the demand error 

23 



T 

Q02153US 

Usl, the operating state of the engine and the like. Even when the 
purification capability of the catalyst converter deviates from the optimal 
state shown by the window 27, the allowable range has a sufficient width to 
allow the catalyst converter to quickly return to the optimal state while 
5 suppressing variations in combustion conditions that may be caused by 
variations in the air-fuel ratio. Therefore, the purification rate of the 
catalyst converter can be kept at high level so that deleterious substances 
of exhaust gas are reduced. 

More specifically, the allowable range is variably updated in 

10 accordance with the determined demand error Usl. For example, the 
allowable range is extended in accordance with deviation of the demand 
error Usl from the allowable range. On the other hand, when the demand 
error Usl is within the allowable range, the allowable range is reduced. 
Thus, the allowable range suitable for the demand error Usl, which defines 

15 the air-fuel ratio necessary to cause the output of the 02 sensor 17 to 
converge to the target value, is established. 

Furthermore, the allowable range is established to be narrower as 
the degree of instability of the output of the 02 sensor 17 becomes higher. 
The allowable range may be established in accordance with the operating 

20 state of the engine including an engine start, an idling state, and 
cancellation of a fuel cut. 

The determined air-fuel ratio error kcmd is added to the base value 
FLAF/BASE to determine the target air-fuel ratio KCMD. The target 
air-fuel ratio KCMD is given to the exhaust system 19 (that is the object to 

25 be controlled), thereby causing the sensor output error Vo2 to converge to 
zero. 

Alternatively, the base value FLAF/BASE of the air- fuel ratio may 
be updated in accordance with the adaptive law input Uadp determined by 
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the sliding mode controller 34 after the completion of the limiting process 
by the limiter 35. More specifically, the base value FLAF/BASE is 
initialized to the stoichiometric air-fuel ratio. If the adaptive law input 
Uadp exceeds a predetermined upper limit value, the basee value 
FLAF/BASE is increased by a predetermined amount. If the adaptive law 
input Uadp is below a predetermined lower limit value, the base value 
FLAF/BASE is decreased by a predetermined amount. If the adaptive law 
input Uadp is between the upper and lower limit values, the base value 
FLAF/BASE is not updated. The base value FLAF/BASE thus updated is 
used in the next control cycle. Thus, the base value FLAF/BASE is adjusted 
to be a central value for the target air-fuel ratio KCMD. 

By performing the above updating process of the base value 
FLAF/BASE in combination with the limiting process, the allowable range 
of the demand error Usl is balanced between positive and negative values. 
It is preferable that the updating process for the base value FLAF/BASE is 
performed when it is determined that the output V62/OUT of the 02 sensor 
substantially converges to the target value Vo2/TARGET and that the 
sliding mode control is in a stable state. 

First decimation filter 

Figure 8 is a block diagram of the first decimation filter 36. A first 
oversampler 51 oversamples the 02 sensor output Vo2/OUT in a shorter 
cycle "n" than the control cycle "k" that is used for calculating the 
manipulated variable Usl for manipulating the air-fuel ratio (that is, "k" is 
the control cycle shown in the above equations). The cycle "n" for the 
oversampling process is, for example, one-fifth of the control cycle "k." The 
oversampled value Vo2_ov is provided to a first low-pass filter 52. 

The first low-pass filter 52 performs a filtering process on the 
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oversampled value Vo2_ov in accordance with the equation (23) to output 
Vo2_ovf. In the equation (23), alovf, a2ovf, a3ovf, bOovf, blovf, b2ovf and 
b3ovf are filtering coefficients that are predetermined with simulation or 
the like. 

5 

Vo2 _ ovf(n) = alovf • Vo2 _ ovf(n - 1) + a2ovf • Vo2 _ ovf(n - 2) + a3ovf • Vout_ ovf(n - 3) 
+ bOovf • Vo2_ov(n) + blovf • Vo2_ov(n - 1) 
+ b2ovf • Vo2 _ ov(n - 2) + b3ovf • Vo2 _ ov(n - 3) 

(23) 

A first downsampler 53 re-samples the filtered value Vo2_ovf in the 
control cycle "k" to output a sample value Vo2_df. 

A method for detecting deterioration of the catalyst based on the 02 
10 sensor output in a certain frequency regions has been proposed. It is 
preferable that the first low-pass filter is designed without invalidating 
such detection of the catalyst deterioration. 

With reference to Figures 9 and 10, a frequency region required to 
detect deterioration of the catalyst will be described. Figure 9 shows a 
15 power spectrum of the 02 sensor output Vo2/OUT (a) when the catalyst is 
new, (b) when the purification rate of the catalyst is sufficient, and (c) when 
the purification rate of the catalyst is insufficient. As seen from Figures 
9(a) through 9(c), the level of the power spectrum of the sensor output 
Vo2/OUT in the frequency region of 3 through 7 Hz varies, which is 
20 indicated by reference number 61. 

Figures 10(a) through 10(c) show a result of filtering the sensor 
output Vo2/OUT shown in Figures 9(a) through 9(c) with a band-pass filter, 
respectively. The power spectrum of the sensor output Vo2/OUT in the 
frequency region of 3 through 7 Hz is emphasized by the filter. As shown by 
25 reference number 62, as the catalyst deteriorates, the power spectrum of 
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the sensor output Vo2/OUT in the frequency regions 3 through 7 Hz 
increases. Thus, by evaluating the sensor output Vo2/OUT in the frequency 
region of 3 through 7 Hz, it can be determined whether the catalyst is in a 
deteriorated state. In order to detect deterioration of the catalyst, the first 
5 low-pass filter 52 is preferably designed not to cut the frequency region of 3 
through 7 Hz. 

Figure 11 shows filter characteristics of the first low-pass filter 52. 
Frequency components necessary to detect deterioration of the catalyst 
exist in a frequency region lower than the line indicated by reference 
10 number 64. The cut-off frequency is set at a frequency sufficiently higher 
than the frequency required for detecting deterioration of the catalyst. 
Thus, the first low-pass filter 52 does not reduce the accuracy of detecting 
the catalyst deterioration. 

15 Effect of the use of the first decimation filter 

Figure 12(a) shows one example of the 02 sensor output 65 sampled 
in the control cycle "k" in accordance with one embodiment of the present 
invention. In the example, the 02 sensor whose response is enhanced was 
used. Figure 12 (a) also shows the output Vo2_df 66 from the first 

20 decimation filter that has the filter characteristics shown in Figure 11. A 
large variation due to chemical noise appears in the 02 sensor output 65, 
as shown in the area 67. Such large variation due to chemical noise is 
removed from the output Vo2_df of the first decimation filter. 

The 02 sensor output 65 shown in Figure 12(b) is the same as that 

25 shown in Figure 12(a). A graph 68 shows Vo2_f obtained by filtering the 
output of the 02 sensor with a low-pass filter that has filter characteristics 
shown in Figure 13. As shown in Figure 13, the low-pass filter has a higher 
cut-off frequency than the frequency required for detecting deterioration of 
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the catalyst. Therefore, the accuracy of detecting deterioration of the 
catalyst is not reduced. It should be noted that the low-pass filter is applied 
to the 02 sensor output that has not been oversampled. 

As shown in the area 69, a large variation occurs in the filtered 
value Vo2_f in accordance with the large variation in the 02 sensor output 
65 caused by chemical noise (although the figure may be hard to see, a 
variation in a convex shape appears in the filtered value Vo2_f in 
accordance with the variation in the 02 sensor output). Thus, in the 
example shown in Figure 12(b), chemical noise included in the 02 sensor 
output cannot be removed. 

The 02 sensor output 65 shown in Figure 12(c) is the same as that 
shown in Figure 12(a). The graph 70 shows Vo2_f obtained by filtering the 
output of the 02 sensor with a low-pass filter that has filter characteristics 
shown in Figure 14. The low-pass filter has a lower cut-off frequency than 
the frequency required for detecting the deterioration of the catalyst, as 
shown in Figure 14. Therefore, the low-pass filter may reduce the accuracy 
of detecting deterioration of the catalyst. The low-pass filter is applied to 
the 02 sensor output that has not been oversampled. 

As shown in the area 71, although a large variation occurs in the 02 
sensor output 65 due to chemical noise, there is little variation in the 
filtered value Vo2_f . The filtered value Vo2_f has a large phase delay 
relative to the 02 sensor output, as clearly seen in the area 72. In the 
example shown in Figure 12(c), not only the accuracy of detecting 
deterioration of the catalyst deteriorates, but also the effect of the improved 
response of the 02 sensor is invalidated. 

Thus, the first decimation filter can remove chemical noise that 
appears in the sensor output Vo2/OUT without causing a phase delay. 

Figure 15 shows one example of the exhaust gas sensor output 
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Vo2/OUT and the target air-fuel ratio error "kernel" in accordance with 
conventional air-fuel ratio control. In the example, the 02 sensor whose 
response is enhanced was used. As shown by reference number 75, 
high-frequency chemical noise appears in the sensor output, which is 
caused by the improved response of the 02 sensor. As shown by reference 
number 76, such chemical noise causes a large variation in the target 
air-fuel ratio error kemd. Such a large variation in the target air-fuel ratio 
error kemd may reduce the purification rate of the catalyst because the 
actual air-fuel ratio of the exhaust gas flowing into the catalyst varies 
according to the target air-fuel ratio error kemd. 

In Figure 16, reference number 81 shows the exhaust gas sensor 
output V02/OUT, reference number 82 shows the first decimation filter 
output V02_df, and reference number 83 shows the target air-fuel ratio 
error kemd, in accordance with one embodiment of the present invention. In 
the example, the 02 sensor whose response is enhanced was used. 
High-frequency chemical noise appears in the sensor output V02/OUT. Such 
chemical noise is removed as shown in the first decimation filter output 
Vo2_df. The use of the first decimation filter output Vo2_df makes the 
target air-fuel ratio error kemd stable. 

Adaptive air-fuel ratio control in accordance with a second embodiment 

Figure 17 shows a block diagram of adaptive air-fuel ratio control in 
accordance with a second embodiment of the present invention. The second 
embodiment differs from the first embodiment in that the first decimation 
filter 36 is removed and in that a second decimation filter 37 is provided. 

The output Vo2/OUT of the 02 sensor 17 is compared with a target 
value Vo2/TARGET. An error Vo2 between the sensor output Vo2/OUT and 
the target value Vo2/TARGET is provided to the controller 31. 
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The output KACT of the LAF sensor 16 is delivered to the second 
decimation filter 37. The second decimation filter 37 oversamples the actual 
air-fuel ratio KACT, filters the oversampled value by a low-pass filtering 
process, and then downsamples the filtered value to output a sample value 
KACT_df. The sample value KACT_df is compared with the target value 
FLAF/BASE. An error "kact" between the sample value KACT_df and the 
target value FLAF/BASE is provided to the controller 31. 

An exhaust system 19, which is an object to be controlled, can be 
modeled as shown by the equation (24), in which V62/OUT is defined as a 
control output and the output KACT of the LAF sensor is defined as a 
control input. 

Vo2(k + 1) = al • Vo2(k) + a2 ■ Vo2(k - 1) + bl ■ kacf(k - d3) 
where Vo2(k) = Vo2 / OUT(k) - Vo2 / TARGET 
kact'(k) = KACT(k) - FLAF/BASE 

(24) 

The sensor output error Vo2 indicates an error between the 02 
15 sensor output Vo2/OUT and the target value Vo2/TARGET. An actual 
air-fuel ratio error kact' indicates an error between the LAF sensor output 
KACT and the base value FLAF/BASE. 

"d3" indicates a dead time in the exhaust system 19. al, a2 and bl 
are model parameters generated by the above-described identifier. 
20 Relation between the LAF sensor output KACT and the second 

decimation filter output KACT_df is expressed as shown by the equation 
(25). 

KACT(k) = KACT _ df (k + d6) (25) 

"d6" indicates a dead time in the second decimation filter 37. The 
25 dead time d6 indicates a time required for the LAF sensor output KACT to 
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be oversampled, filtered by the low-pass filtering process, and then 
dwonsampled. The dead time d6 is, for example, one control cycle (that is, 
d6=l). 

A system comprising the exhaust system 19 and the second 
decimation filter 37 can be obtained from the equations (24) and (25), as 
shown by the equation (26). 

Vo2(k + 1) = al - Vo2(k) + a2 • Vo2(k - 1) + bl • kact(k - d3 + d6) 
= al - Vo2(k) + a2 . Vo2(k - 1) + bl • kact(k - dl) 
where Vo2(k) = Vo2 / OUT(k) - Vo2 / TARGET 
kact(k) = KACT _ df (k) - FLAF/ BASE 
dl = d3 - d6 

(26) 

Thus, incorporation of the second decimation filter 37 decreases the 
dead time in the exhaust system. 

On the other hand, the air-fuel ratio manipulating system 18 
extending from the ECU 5 through the engine 1 to the LAF sensor 16 can 
be modeled as shown by the equation (27). 

kact'(k) = kcmd(k-d4) 

where kact'(k) = KACT(k) - FLAF /BASE 
kcmd(k) = KCMD(k) - FLAF /BASE 

(27) 

The target air-fuel ratio error "kcmd" indicates an error between the 
target air-fuel ratio KCDM and the base value FLAF/BASE. "d4" indicates 
a dead time in the air-fuel ratio manipulating system 18. 

A system comprising the air-fuel ratio manipulating system 18 and 
the second decimation filter 37 can be obtained from the equations (25) and 
(27), as shown by the equation (28). 
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kact(k) = kcmd(k - d4 - d6) 
= kcmd(k-d2) 

where kact(k) = KACT _ df (k) - FLAF / BASE 
kcmd = KCMD - FLAF /BASE 
d2 = d4 + d6 

(28) 

Thus, incorporation of the second decimation filter 37 increases the 
dead time in the air-fuel ratio manipulating system. 

The air-fuel ratio manipulating system 18 may be included in the 
object to be controlled by the adaptive air-fuel ratio control. In this case, 
the model equation for the controlled object is expressed based on the 
equations (26) and (28), as shown by the equation (29). The dead time "d" is 
a total dead time in a system comprising the air-fuel ratio manipulating 
system 18, the exhaust system 19, and the second decimation filter 37. 

Vo2(k + 1) = al - Vo2(k) + a2 - Vo2(k - 1) + bl • kcmd(k - dl - d2) 
= al - Vo2(k) + a2 ■ Vo2(k - 1) + bl • kcmd(k - d) 
where d = dl + d2 = d3 - d6 + d4 + d6 = d3 + d4 

(29) 

Thus, the dead time in the system comprising the air-fuel ratio 
manipulating system 18, the exhaust system 19, and the second decimation 
filter 37 is equal to a sum of the dead time d3 in the exhaust system 19 and 
the dead time d4 in the air-fuel ratio manipulating system 18. The dead 
time d6 in the second decimation filter has no influence on the system. 

The structure of the controller 31 in the second embodiment is the 
same as that in the first embodiment. The dead time dl, d2 and d shown in 
the above equations regarding the identifier 32, the estimator 33 and the 
sliding mode controller 34 of the controller 31 are (d3~d6), (d4+d6) and 
(d3+d4), respectively. 
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Second decimation filter 

Figure 18 is a block diagram of the second decimation filter 37. A 
second oversampler 55 oversamples the LAF sensor output KACT in a 
shorter cycle "n" than the control cycle k that is used for calculating the 
5 manipulated variable Usl for the air-fuel ratio control (that is, "k" is the 
control cycle used in the above equations). The cycle "n" for the 
oversampling process is, for example, one-fifth of the control cycle "k." The 
oversampled value KACT_ov is provided to a second low-pass filter 56. 

The second low-pass filter 56 filters the oversampled value 
10 KACT_ov in accordance with the equation (30) to output KACT_ovf. In the 
equation (30), alovf , a2ovf , a3ovf , bOovf , blovf , b2ovf and b3ovf are filter 
coefficients predetermined with simulation or the like. Filter 
characteristics of the second low-pass filter 56 is shown in Figure 19. 

KACT_ ovfln) = alovf KACT_ ovf(n - 1) + a2ovf-KACT_ ovfln - 2) + a3ovf-KACT_ ovfln - 3) 
+ bOovf -KACT_o^n) + blovf -KACT_ov(n - 1) 
15 +b2ovf-KACT_o\<n-2) + b3ovf.KACT_o\{n-3) 

(30) 

A second downsampler 57 re-samples the filtered value KACT_ovf in 
the control cycle "k" to output a sample value KACT_df. 

20 Effect of the use of the second decimation filter 

In Figure 20, reference number 91 shows the actual air-fuel ratio 
error kact detected by the LAF sensor, and reference number 92 shows the 
target air-fuel ratio error kcmd calculated based on the error kact in 
accordance with conventional air-fuel control. In the example, the LAF 
25 sensor whose detection range is expanded was used. Reference number 93 
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shows the model parameter bl calculated by the identifier 32. As described 
above, the model parameter bl acts to identify correlation between the 
air-fuel ratio error kact, which is an input of the modeled control system, 
and the sensor output error V62, which is an output of the modeled control 
system. 

As seen, for example, in the area 94, even if the actual air-fuel ratio 
error kact continuously varies within a range between -P2 and +P2, the 
actual air-fuel ratio error kact takes one of the values of -P2, -PI, 0, +P1 
and +P2 due to the low resolution of the LAF sensor. For example, an 
actual air-fuel ratio error that has an intermediate value between 0 and 
+P1 cannot be detected. The actual air-fuel ratio error kact is often detected 
to be continuously zero as indicated by a flat portion shown in the area 94. 

If such a flat portion is detected continuously, the identifier 
determines that there is no correlation between the air-fuel ratio error kact 
and the sensor output error Vo2. If it is determined that there is no 
correlation, the identifier decreases the model parameter bl as shown by 
reference number 93. Thus, the identification accuracy of the model 
parameter bl deteriorates. When the value of the model parameter bl 
becomes small, the sliding model control reaches an over-corrected state. 
This state increases variation in the target air-fuel ratio error kcmd, as 
indicated by the region 95 of the graph 92, thereby reducing the 
purification rate of the catalyst. 

In Figure 21, reference number 97 shows the actual air-fuel ratio 
error kact detected by the LAF sensor, and reference number 98 shows the 
output KACT_df from the second decimation filter 37 in accordance with 
the second embodiment of the present invention. The LAF sensor whose 
detection range is expanded was used. It is seen that values that cannot be 
detected by the LAF sensor (that is, values below the resolution limit of the 
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LAF sensor) are estimated by applying the second decimation filter to the 
actual air-fuel ratio KACT. There is no continuous flat portion in KACT_df. 
Therefore, the identifier determines that there is correlation between the 
actual air-fuel ratio error kact and the sensor output error Vo2. The value 
5 of the identification parameter bl becomes stable as shown by a graph 99. 

Thus, the second decimation filter compensates the shortage of 
resolution of the LAF sensor. The cycle "n" used for the oversampling 
process in the second decimation filter 37 may be the same as the cycle 
used for the oversampling process in the first decimation filter 36. 
10 Alternatively, the cycle "n" used in the second decimation filter 37 may be 
different from the cycle used in the first decimation filter 36. 



Adaptive air-fuel ratio control in accordance with a third embodiment 

Figure 22 shows a block diagram of adaptive air-fuel ratio control in 
15 accordance with a third embodiment of the present invention. In the third 

embodiment, both of the first decimation filter 36 and the second 

decimation filter 37 are provided. 

The output Vo2/OUT of the 02 sensor 17 is delivered to the first 

decimation filter 36. The output Vo2_df from the first decimation filter 36 is 
20 compared with the target value Vo2/TARGET. An error Vo2 between the 

sample value Vo2_df and the target value Vo2/TARGET is supplied to the 

controller 31. 

The output KACT of the LAF sensor 16 is delivered to the second 
decimation filter 37. The output KACT_df from the second decimation filter 
25 37 is compared with the target value FLAF/BASE. An error kact between 
the sample value KACT_df and the target value FLAF/BASE is supplied to 
the controller 31. 

The exhaust system 19, which is an object to be controlled, can be 
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modeled as shown in the equation (31), in which Vo2/OUT is defined as a 
control output and the LAF sensor output KACT is defined as a control 
input. 

Vo2'(k + 1) = al - Vo2'(k) + a2 • Vo2'(k - 1) + bl • kact'(k - d3) 
where Vo2' (k) = Vo2 / OUT(k) - Vo2 / TARGET 
kact 1 (k) = KACT (k) - FLAP / BASE 

(31) 

A sensor output error Vo2' indicates an error between the 02 sensor 
output Vo2/OUT and the target value Vo2/TARGET. An actual air-fuel ratio 
error kact' indicates an error between the LAF sensor output KACT and the 
base value FLAF/BASE. 
10 "d3" indicates a dead time in the exhaust system 19. al, a2 and bl 

are model parameters generated by the above-described identifier. 

Relation between the 02 sensor output Vo2/OUT and the output 
Vo2_df of the first decimation filter is represented by the above equation (2). 
Relation between the LAF sensor output KACT and the output KACT_df of 
15 the second decimation filter is represented by the above equation (25). 

A system comprising the exhaust system 19, and the first and the 
second decimation filters 36 and 37 can be obtained based on the equations 
(2), (25) and (31), as shown by the equation (32). 



Vo2(k + 1) = al • Vo2(k) + a2 • Vo2(k - 1) + bl ■ kact(k - d3 - d5 + d6) 
= al • Vo2(k) + a2 - Vo2(k - 1) + bl • kact(k - dl) 
where Vo2(k) = Vo2 _ df (k) - Vo2 / TARGET 
kact(k) = KACT _ df (k) - FLAF / BASE 
dl = d3 + d5-d6 

(32) 

If a dead time d6 in the second decimation filter 37 is equal to a 
dead time d5 in the first decimation filter 36, a dead time of the exhaust 
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system in which the first and the second decimation filters 36 and 37 are 
provided is d3. That is, incorporation of the first and the second decimation 
filters 36 and 37 has no influence on the exhaust system. 

On the other hand, the air-fuel ratio manipulating system 18 
extending from the ECU 5 through the engine 1 to the LAF sensor 16 is 
modeled as shown by the above equation (27). A system comprising the 
air-fuel ratio manipulating system 18, the first and the second decimation 
filters 36 and 37 is expressed by the above equation (28). The dead time in 
the first decimation filter has no influence on the air-fuel ratio 
manipulating system 18. 

The air-fuel ratio manipulating system 18 may be included in the 
object to be controlled by the adaptive air-fuel ratio control. In this case, 
the model equation of the object to be controlled is expressed based on the 
equations (32) and (28), as shown by the equation (33). A dead time d is a 
total dead time in a system comprising the air-fuel ratio manipulating 
system 18, the exhaust system 19, and the first and second decimation 
filters 36 and 37. 

Vo2(k + 1) = al • Vo2(k) + a2 ■ Vo2(k - 1) + bl • kcmd(k - dl - d2) 
= al • Vo2(k) + a2 ■ Vo2(k - 1) + bl • kcmd(k - d) 
d = dl + d2 = d3 + d5 - d6 + d4 + d6 = d3 + d4 + d5 

(33) 

Thus, the dead time in the system comprising the air-fuel ratio 
manipulating system 18, the exhaust system 19 and the first and second 
decimation filters 36 and 37 is equal to a sum of the dead time d3 in the 
exhaust system 19, the dead time d4 in the air-fuel ratio manipulating 
system 18, and the dead time d5 in the first decimation filter. The dead 
time d6 in the second decimation filter has no influence on the system. 

The structure of the controller 31 in the third embodiment is the 
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same as that in the first embodiment. The dead time dl, d2 and d shown in 
the above equations regarding the identifier 32, the estimator 33 and the 
sliding mode controller 34 of the controller 31 are (d3+d5*d6), (d4+d6) and 
(d3+d4+d5), respectively. 

According to the third embodiment, chemical noise is removed from 
the output of the 02 sensor. The shortage of resolution of the LAF sensor 
can be compensated. 

Other embodiment of the controller 

In the above embodiments, the controller 31 (Figures 4, 17 and 22) 
may perform other response assignment control instead of the sliding mode 
control. 

Figure 23 shows another embodiment of the controller. The 
difference from Figure 4 is that the controller 31 comprises a DSM 
controller 38 for performing a AX (delta-sigma) modulation algorithm 
instead of the sliding mode controller. 

Figure 24 shows a block diagram of the DSM controller 38. Since the 
AX modulation algorithm is conventionally used, detailed description is 
omitted. An amplifier 101 multiplies a sensor output error Vo2 , which is 
estimated by the estimator 33, by a gain "-G" to output a reference signal 
"r." A subtracter 102 calculates a difference between the reference signal r 
and the AX modulation signal u"(k-l) calculated in the previous cycle to 
output a differential signal 8(k). 

An integrator 103 adds the differential signal 8(k) calculated in the 
current cycle to the integrated signal a f (k-l) calculated in the previous cycle 
to output an integrated signal o'(k) for the current cycle. The sign function 
104 determines +/- of the integrated signal a'(k) to output a signal u"(k). An 
amplifier 105 multiplies the signal u M (k) by a gain "F" to output a AX 
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modulation signal u(k). The A£ modulation signal u(k) is a signal 
indicating the target air-fuel ratio error kcmd. The target air-fuel ratio 
KCMD is calculated by adding the base value FLAF/BASE to the AX 
modulation signal u(k). 
5 Alternatively, the DSM controller 38 may be configured to perform a 

A (delta) modulation algorithm, which does not contain the integrator. The 
DSM controller 38 may also be configured to perform a £A (sigma-delta) 
modulation algorithm, in which the integrator is followed by the subtractor. 

In yet another embodiment, the controller 31 comprises both of the 

10 sliding mode controller 34 and the DSM controller 38. The controller 31 can 
switch between sliding mode control (or another response assignment 
control) and AX modulation control (or £A modulation control, A modulation 
control) in accordance with the operating state of the engine. For example, 
the air-fuel ratio is controlled by the A£ modulation controller when the 

15 load of the engine is low, and the air-fuel ratio is controlled by the sliding 
mode controller when the load of the engine is high. 

Control flow 

Figure 25 shows a flowchart of the adaptive air-fuel ratio control in 
20 accordance with the above third embodiment of the present invention. The 
cycle "n" for the oversampling process in the first decimation filter is the 
same as the cycle for the oversampling process in the second decimation 
filter. This routine is performed in cycle "n." 

In step Sill, the output Vo2/OUT of the 02 sensor is sampled and 
25 then filtered by the low-pass filter, to determine Vo2_ovf. In step S112, the 
LAF sensor output KACT is sampled and then filtered by the low-pass filter 
to determine KACT_ovf. As described above, since the cycle "n" is shorter 
than the cycle "k" that is used for calculating the manipulated variable Usl 
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for the air-fuel ratio control, the 02 sensor output Vo2/OUT and the LAF 
sensor output KACT are oversampled. Steps Sill and S112 may be 
performed in parallel. 

In step S113, it is determined whether a counter CPRISM is zero. 
The initial value of the counter CPRISM is "k/n-l." For example, when "n" 
is ten milliseconds and "k" is fifty milliseconds, the initial value of the 
counter CPRISM is "4." When the value of the counter CPRISM is not zero, 
the counter is decremented by one in step S114, exiting the routine. If the 
counter CPRISM is zero when the routine is reentered, the process 
proceeds to step S115, in which the counter is reset. Thus, steps Sill and 
S112 are carried out in cycle "n" while steps S115 through S122 are carried 
out in cycle "k." 

In step S116, it is determined whether a permission flag is one. The 
permission flag is a flag that is to be set to one when the execution of the 
air-fuel ratio control is permitted. For example, the execution of the air-fuel 
ratio control is not permitted when lean operation is being performed in the 
engine or when operation for retarding the ignition timing is being 
performed. When the permission flag is zero, the target air-fuel ratio error 
kcmd is set to a predetermined value in step S123, and then the process 
proceeds to step S124. 

When the permission flag is one, Vo2_ovf filtered in step Sill is 
sampled in cycle "k" to determine Vo2_df. KACT_ovf filtered in step S112 is 
sampled in cycle "k" to determine KACT_df. 

In step S118, the calculation process by the identifier is performed 
to determine the model parameters al, a2 and bl as described above. In 
step S119, the calculation process by the estimator is performed to 
determine the estimated value Vo2 as described above. In step S120, the 
manipulated variable Usl is calculated as described above. 
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In step S121, the stability of the sliding mode control is determined. 
For example, it is determined that the sliding mode control is stable when a 
stability determination base parameter Pstb (=a(k + d)-Aa) is equal to or 
less than zero. Pstb is equivalent to the time-differential of the Lyapunov 
function a 2 /2 concerning the switching function a . The state in which the 
function value Pstb is equal to or less than zero indicates a state in which 
the value of the switching function a converges to zero or is converging to 
zero. The state in which the function value Pstb is greater than zero 

indicates a state in which the value of the switching function a is leaving 
zero. 

In step S122, the limitation process is performed by the limiter to 
determine the target air-fuel ratio error kcmd. In step S124, the base value 
FLAF/BASE is added to the target air-fuel ratio error kcmd to determine 
the target air-fuel ratio KCMD. If the base value FLAF/BASE is established 
as described above, a step for establishing the base value may be performed 
after step S124. 

The invention may be applied to an engine to be used in a 
vessel-propelling machine such as an outboard motor in which a crankshaft 
is disposed in the perpendicular direction. 
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